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Abstract 

In this work we intend to study a class of time-dependent quantum systems with non- 
Hermitian Hamiltonians, particularly those whose Hermitian counterpart are important 
for the comprehension of posed problems in quantum optics and quantum chemistry, 
which consists of an oscillator with time-dependent mass and frequency under the action 
of a time-dependent imaginary potential. The propagator for a general time-dependence 
of the parameters and the wave-functions are obtained explicitly for constant frequency 
and mass and a linear time-dependence in the potential. The wave-functions are used to 
obtain the expectation value of the Hamiltonian. Although it is neither Hermitean nor 
PT symmetric, the case under study exhibits real values of energy. 
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In the last few years, the so called PT symmetric systems introduced in the seminal paper 
of Bender and Boettcher [1] has attracted much attention. They consist of non-Hermitian 
Hamiltonians with real eigenvalues, which however exhibit parity and time-reversal symmetries 
or , in other words, when one makes the space-time inversion (P : x — > — x; T : t — > —t, i — > 
— i). In fact, following reference [2], the time-reversal operation Ti/) {x,t) = tp* (x, — t), keeps 
the Schroedinger equation covariant even for time-dependent Hamiltonians [3], [4]. Besides, 
there exist other classes of Hamiltonians having real spectra without being PT symmetric, 
as can be seen, for instance, in references [5]- [6], and systems where the PT symmetry is 
spontaneously broken [7], leading to complex energy eigenvalues. 

Most of the papers dedicated to the subject refers to time-independent systems and station- 
ary states. Despite an enormous list of works devoted to the developing and the understanding 
of this new kind of situation [1, 5-12], as far as we know, no one has discussed the case of 
time-dependent systems, which are very important for the comprehension of many problems in 
quantum optics, quantum chemistry and others areas of physics [13]- [28]. In particular we can 
cite the case of the electromagnetic field intensities in a Fabry-Perot cavity [16]. In fact this 
kind of problem constitutes a line of investigation which still attracts the interest of physicists 
[25]-[28]. 

Here we intend to solve a class of time-dependent quantum system with non-Hermitian 
Hamiltonians, particularly those which consist of an oscillator with time-dependent mass and 
frequency, plus a linear non-Hermitian term with time-dependent coupling. After that, we 
will discuss the effect of that time-dependence on the reality of the spectrum, versus the PT 
symmetry, by treating a particular case. The system to be studied is represented by the one- 
dimensional non-Hermitian Hamiltonian, namely 

which in principle can preserve the PT symmetry, depending on the time-dependence of their 
parameters. In fact it is easy to realize that the PT symmetry imposes that A {t) must be an 
even function of time for m (t) and cu^ (t) symmetric under time-reversal. 
The corresponding Schroedinger equation is written as 
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= (2) 



where m (t), ui it) and A it) are time-dependent functions (We sei h — 1 for a while). 

Usually the time-independent forced harmonic oscillator is solved by means of a translation 
of the coordinates. When we deal with a time-dependent forced harmonic oscillator one can 
check that a trivial translation of the coordinates x = y — i 2m{t)% (t) "^oes not lead to a harmonic 
oscillator with a time-dependent mass and frequency because we are still left with a first order 
derivative term in y and an additional purely time-dependent term (see Eq. (5) below). 

To show that the Schroedinger equation (2) can be mapped into an evident exactly solvable 
equation we appeal to a generalization of transformations, which envolves a change of variable 
and time reparametrization, used before in several time-dependent systems [17], [19], [21], 
namely 

X = s{t) y + ir]{T) , (3) 
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where r is a single-valued function related to the original time t by 

r{t) = 



dt 



(4) 



In the transformation (3), ?7 is a real function and y is necessarily complex to keep the original 
physical variable x real. After the transformation of variables we are led to the following 
equation 
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At this point, in order to get the Schroedinger equation in the new variables, we introduce a 
wavefunction redefinition 

ij{y,T) = exp[if {y,T)]a{y,T) , (6) 
which once substituted in the above equation gives us 
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where the prime denotes differentiation with respect to r. Now we choose conveniently the 
arbitrary function f{y,T) to guarantee that the coefficient of the term |^ vanishes identically. 
In doing so we get 

df ^ , ■ ^ 

— = m u, s s y + tmsuri, (8) 
dy 

which has as a general solution 



f{y,^)= ^muss y^ +imsiir}y+ (r) , 



(9) 



with (r) an arbitrary function of the rescaled time r, which is still to be determined by some 
convenient imposition. Now, substituting this function into the equation (7), one can rewrite 
it as 



' d 1 

^A* — + 



-I fj, 



dr ' 2m dy"^ 2 
X 



muj s d ,2 

\- — ymiis s) — iim s 

jjL dr 



y'+ 



mtu's 



r] + 



muj' 



dr 



is dfr I f 2X ^ 

+ o/^ + r][r]+ 

2 s dr 2 \ mcu^ 



1 '1 
miJ, r) 



a{y,T) = 0. 



We choose the (r) and impose some constraints on the functions s(t), ?7(r) and //(r) in order 
to leave the effective potential depending only on y^. That is fr (r) must to be such that 
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whose solution can be cast into the form below without loss of generality 
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Besides s{t), r}{T) and //(r) have to satisfy the equation 
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which allow us to write the corresponding Schroedinger equation in the form 
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Until now, the functions s (r) and (r) remains arbitrary since cqTiation (12) constrains ?7(r) in 
terms of s (r) and /i (r). Then we use this arbitrariness to get a simpler equation to this problem. 
One obvious possibility is choosing them in order to get a time-independent Schroedinger 
equation for a {y, r) 



f.d 1 a' 
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(15) 



where it can be seen that the time dependence was factorized, implying that 



ms^jjL = mo = const.; fa;^ + fl^) = mo<^n = const. 



(16) 



Prom the above equations one can see that we have transformed the original problem into 
that of an usual harmonic oscillator with constant mass and frequency. At this point we are 
able to calculate the propagator for the system, once it is a special solution of the Schroedinger 
equation, subject to the restriction 



lim K (x, t; Xq, to) — 5 {x — Xq) . 

Now, remembering that in this case the wavefunction is written as 

i;(x,t) = e'f^y'^^a(y,T) 



(17) 



(18) 



4 



and finally using that 

/oo 
K {x,t;xo,to) 'ilj{xo,to) , 
-oo 

it is easy to conclude that the propagator will come from [19] 



K{x,t;xo,to)^ e'f^y'^^K{x,T;xo,To) e-^^*(^°'^° 



y={x -i »7(T))/s(r), r=r(i) 



(19) 



(20) 



where f* {yo,To) stands for the complex conjugate evaluated at initial time and coordinate. 
Now using the well known expression of the propagator of the harmonic oscillator [29], one 
obtains 



K (x,t;xo,to) = 
mo u!o 



+ 
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, (21) 

sm [Uo [T-To)\ ) y^^^ _ . ^(^t))/s(T), T=T(t) 

where it can be observed that we recall the dependence on h, and use the identification /is' = s, 
where the dot denotes differentiation with respect to t. Furthermore one has a nonformal 
solution provided that an exact solution of (12) and (16) is known [17] -[2 5], [2 8]. 

Let us now obtain an expression for the wavefunctions directly from the above propagator. 
For this we make use of the Mehler's formula [30] 



exp 
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with a = ^J ^°^° yo, b — ..J^^y, c — exp {—i [u>o (r — Tq)]}. Moreover, the propagator can be 
written in its spectral decomposition form 



K {y, r; yo, tq) = XI C {y, ^ {yo, tq) , 



(23) 



n=0 



and from which the wavefunctions can be devised as 

1 
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(24) 



y={x -i r){t))/s(T), T=T{t) 

In the above expressions i/„ stands for the Hermite polynomial. 

Let us now treat a simpler particular case, which can be used to get a deep understanding 
of the problem without unnecessary mathematical complications namely, the case of a usual 
harmonic oscillator with a time-dependent PT- violating linear potential, 

d'^ijj{x,t) 



2 m dx"^ 



+ 



m uj^ x^ + iatx 



ip {x,t) — ih 



dilj {x, t) 
dt ^ 



(25) 
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where a is an arbitrary real constant. The coordinate and time transformations given in equa- 
tions (3) and (4) can be simphfied by choosing s(r) = yu(r) = 1, (r = t) and rj{t) = — 
Then one ends up with the Schroedinger equation (15) where coq — uj and mo = m, whose 
solutions in the transformed coordinate looks like simply as 



At this point wc are able to calculate the energy expectation value of the system. For this 
we calculate, as usual, the expression 

^ ' J^^dxi;*{x,t)i;{x,t) • ^ ^ 

After straigthforward calculations one can get 



where 



n + -j fiuj —h 7(t) — muj f] 
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We do not have a general expression for the energy expectation value for an arbitrary 
eigenstate, but we have verified that they are always real. In fact, the calculation for a particular 
state can be done straightforwardly and we have, for instance 



(o|£;|o) = -Tiuj 

{1\E\1) = ^na; 

This way we have proven that the reality of the time-dependent energy is mantained despite 
of the non-Hermiticity of the Hamiltonian. For a non-Hermitian and PT- violating Hamiltonian 
this is a very intriguing conclusion and demands further investigations on the origin of the 
reality of the energy in such systems, which can also be strongly considered as candidates 
for realizations of non-Hermitian interactions. Furthermore, one can note that even with the 
mapping into a Hermitian time-independent system, Eq. (15), it is not trivial to verify that the 
original time-dependent system has real energies. In fact, it is necessary to use the wavefunction 
(6) to calculate the energy expectation value. We are still studying cases where the mass and 
frequency are constants, looking for the more general time-dependent forces, including PT 
symmetric ones, that present real spectrum. 
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